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Abstract 

Using a generalized proper-time method, we obtain expressions for the fermion 
density and the QED effective Lagrangian for an external magnetic field at 
finite chemical potential. The effective Lagrangian and the density are here 
written in terms of elementary functions, summed over a finite number of filled 
Landau levels. 
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The study of finite temperature and density quantum electrodynamics (QED) with a 
nonvanishing average magnetic field is of considerable interest as it may be associated with 
the electron-positron plasma in compact stellar objects (e.g. neutron stars and magnetic 
white dwarfs), where the fermion density and the magnitude of the magnetic field may 
be extremely high (see e.g. Refs. |l], U). 

The QED thermodynamical potential at finite temperature and density with a static 
uniform magnetic field was calculated already 25 years ago in Ref . @ , and using a gener- 
alization of Fock-Schwinger's proper-time method for T,/i ^ (where /i is the chemical 
potential) later in Ref. [|J. The interest to this problem was renewed in Ref. [||, where 
an elegant generalization of Fock-Schwinger's proper-time method in the case of nonzero 
magnetic field and chemical potential was made. Using a real-time thermal formalism, 
the results of Ref. || were completed and generalized to finite temperature in Ref. ||. 
The expressions obtained for the effective action in the above cited references were rather 
complicated and did include some proper-time like integrals and/or infinite sums. 

Here we want to demonstrate that using the formalism, elaborated in Ref. |J we 
may move forward and for the finite density QED with an external magnetic field at 
zero temperature show that it is possible to obtain simple expressions for the fermion 
density and the effective Lagrangian. The effective Lagrangian is here written in terms 
of elementary functions as a sum over a finite number of (partially) filled Landau levels, 
and agrees with the zero temperature limit of the fermion partition function. As an 
application we confirm that the magnetisation obtained from this effective action does 
exhibit the relativistic de Haas-van Alphen effect |3|, || [7|] . 

We shall here consider finite density QED with a nonvanishing average magnetic field. 
Including the chemical potential the corresponding Lagrangian reads 

£ = -^F^ + ^-e4- 7o //-m)V , (1) 

where we have chosen the gauge = (0, —X2B, 0, 0). 

In order to calculate the one-loop correction to the effective Lagrangian, / <i1r£ cfr = 
— zln Det(i^ — eJji. — 70/1 — m), we shall first evaluate the fermion density p = We 
may then reconstruct the effective Lagrangian according to 

£ eS (B,ri=£ eS (B)+fr s (B,n) , (2) 

where 

£ eS (B,n)= f p{B,n')dn' , (3) 
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is the contribution due to the finite density, and 
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is the well known vacuum part of the effective Lagrangian in the purely magnetic case 
We may rewrite the expression for the fermion density as 



p = i tr (7oG(:r = x')) = i 
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where the trace is over spinor indices only, and G(x, x') is the fermion Green's function 
in configuration space. We shall use the expression obtained in Ref. || for the Green's 
function in momentum space 
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—i9((po + /x)signp ) J ds exp < is 
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{(1 - 7i7 2 tan(eSs)) 7 3 p|| - 70 (p + p) - m + {1 + tan 2 (e.Bs)}(7ipi + 72P2)} (6) 

where p» and p± are the modulus of the momenta parallel and perpendicular to the 
magnetic field, respectively. A different ie-prescription here arise for \fi\ > \m\ as the 
rules for passing poles in the fermion Green's function are changing H [| , since the Dirac 
sea is filled up to the energy //. 

Using Eq.(^j) in Eq.(^), the identity 9{x) = 1 — 9(—x), and performing a trivial change 
of variables, we get 



p( B ,p) = J <&P J o dsp exp 
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where the first term on the right-hand side is vanishing. 
Performing the momentum integration in Eq.(^), we get 



p(B,fi) = 23?e 



e 3ivr/4 

8^ 



ds 



(eBs) cot(eBs) exp{is(/i 2 — m 2 + is)} 



that may be obtained from the expression for the effective Lagrangian in Ref. ||, keeping 
the is prescription. This is prescription tells us that the proper-time integral actually 
is to be performed slightly below the real axis. When closing the contour to obtain an 
exponentially decreasing integrand instead of an oscillating one, the poles of cot(eBs) 
will thus be encircled when the contour is closed in the upper half-plane for fi 2 > m 2 . 
The sum over the residues at these poles will form the vanishing temperature limit of the 
"oscillating" part of the effective Lagrangian, in agreement with Ref. ||. 

The proper-time integral in Eq.@, or the corresponding integral after the above- 
described Wick rotation, cannot be performed analytically. Instead, we shall here perform 
the proper-time integral in Eq. (|7]) before the integral over the momentum. This equation 



may be rewritten as [1C 
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where we have defined 

Hp) = J ds exp < is 
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For m 2 > p§ — p| we may close the integration contour in the lower half plane 
I(p) = — i J ds exp < — s 



2 2,2,2 tanh(e5s) 
m ~Po+P\\ +P± l£ 
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The integral in Eq.([TI]) is diverging as s —>■ oo for m 2 < p\ — p 2 . Changing the variable 
of integration to z — tanh(e-Bs) (eB > 0), and defining Q = (m 2 + p 2 — p^ — ie) /2eB, 
we may rewrite Eq.(|TTl) as JlTJ 



-1+0(1 + exp 



(12) 



that has a singularity at z = 1. From the left- and right-hand sides of Eq.([T|) one may 
now extract the first k terms of the Taylor expansion of (l + z)~ 1- exp{—p 2 ± z/ eB} around 
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z = 1 (Cauchy method), 

Hp) + ~bY. dza n (p)( 
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where we have defined 



a n (p) = L_i_ (i + z )-i-««p{-^/ C B} , . (14) 



n\ dz r 
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Performing the trivial integration in the left-hand side of Eq.flToD, under the assumption 
that Q + n > 0, we get 

(15) 

The integral in the right-hand side of Eq . fll5|) is convergent in the halfplane dte(Q) > 
— (k + 1). Taking limit k — > oo we see that Eq . ((TBI) is an analytical continuation of I(p) 
on the whole complex plane p\ — pjj , excluding the points — Py = 777,2 + "2eBn, n = 
0, 1,2, . . ., where J(p) has simple poles (these poles are nothing but the familiar rela- 
tivistic Landau levels). 

Substituting the expression for I(p) rewritten as in Eq. fll5D into Eq.@ we see that 
the regular parts of I(p) and I*(p) cancel, and thus 
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P( B > A*) = 0-4 / Po^Po / ^pE^ 

J n=0 
On(p) <(P) 



(16) 



Po ~ Py — 7772 — 2e -Bn + Po — £>| — m2 — 2eBn — ie 

Using the identity (x ± ie)^ 1 = p(x _1 ) =F in5(x), we see that also the principal values 
cancel, and the only nonvanishing contribution comes from the poles 

p(B,fx) = — Podpo d 3 pJ2a n (p)5(pl- pf\- m 2 - 2eBn) 
vr jo J n=0 
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I roo roo 00 

= T~2 / dp W / rf ^J-) H a n(P(n))0(^ 2 -p\-m 2 - 2eBn) , (17) 

where £>(„) denotes the four- momentum, such that (po(n)) 2 = m2 + G°||(n)) 2 + 2eBn. The 
Heavyside step-functions here does accordingly describe the number of (partially) filled 
Landau levels. For example, for < p 2 — m 2 < 2eB, only the lowest Landau level (n = 0) 
contributes to the density 

I roo roo 

pa{B,n) = — J d P\\ J Q d(p\)a (p ) 6(p 2 - p\ - m 2 - 2eBn) 
pB i 

= ^vV 2 -™ 2 ^ 2 -m 2 ) . (18) 
For < /i 2 — m 2 < 4e£>, we get similarly p(B,p) = p (B,p) + pi(B,p), where 

Pl (B,p) = _ m 2_ 2eB 9{fl 2 _ m 2_ 2eB) j~> d(pl)?± exp L&\ 

eB 
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~—^p 2 -m 2 -2eBQ{p 2 -m 2 -2eB) . (19) 



The general expression for the fermion density in a static uniform magnetic field B, may 
thus be written as a sum over a finite number of occupied Landau levels 



p(B,p) = J2 Pn(B,p) , (20) 



n=0 



where [. . .] denotes the integral part. The contribution from the n-th Landau level is 

eB 



Pn (B,p)=b n —^p 2 -m 2 -2eBn , (21) 
where we have defined 

b n = 2-5 nfi , (22) 

since the lowest Landau level (n = 0), unlike the higher levels, only contains fermions 
with one projection of the spin; as found in the above examples of po and p%, cf. Eq.(p6|). 

We have thus found an expression for the fermion density in a nonvanishing average 
magnetic field, in terms of elementary functions in a discrete finite sum over filled Landau 
levels. This result may be well understood from the index theorem approach fll2| . The 



density (fermion number) depends on the difference of numbers of filled positive and 
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negative energy levels, which in this case is the (semidescrete) number of Landau levels 
in the interval [0,/i]. 

In the limit of vanishing magnetic field, eB — > 0, the Riemann sum of Eq.(|2T|) may be 
rewritten as an integral 

p( /i ) = J- dx( f i 2 -m 2 -x) 1 / 2 6( l 2 2 -m 2 ) = U 2 - m 2 ) 3/2 fl^ 2 - m 2 ) (23) 
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which is the familiar expresion for the fermion density. 

In Figure 1 we show the density as a function of the chemical potential for fixed 
magnetic field, and in Figure 2 the density is given as a function of the magnetic field for 
fixed chemical potential. We see that the density is showing an oscillating behaviour as 
consecutive Landau levels are passing the Fermi level. 

Integrating Eq . (|2~0D with respect to the chemical potential, we find the part of the 
effective Lagrangian due to the finite density as 



2eS 



£ cS (B,f,)= C n{B^) , (24) 

n=0 

where the contribution from the n-th Landau level is 

C n (B, ») = b n ^- 2 LvV 2 - m 2 - 2eBn - (m 2 + 2eBn) In f £ + V f ~ if Z ) . 

4:TT Z [ v V y/m 2 + 2eBn J J 

(25) 

Here we have used the zero temperature proper-time method to calculate the fermion 
density and effective Lagrangian. It is noteworthy to compare with the approach of 
quantum statistical mechanics. By comparing the generating functional for fermionic 



Green's functions in imaginary time |T3[ , with the partition function in the grand canonical 



ensemble (Z), we find that C efC = 4y In Z. The relativistic fermion energy-levels in a static 



uniform magnetic field are found as |T4] 



E k ,x(p\\) = sJm?+ V \ + 2eB{k + \-l) , (26) 

where k = 0, 1, 2, . . . corresponds to the quantized orbital angular momentum, and A = 
1,2 describes the projection of spin. Using the ordinary relativistic dispersion law to 
reintroduce the momenta orhogonal to the magnetic field, we find the density of states 
V eB / (2tt) 2 , and readily obtain 

C cS (B,fi,T) = 



1 eB 

/?(2tt) 2 



OO 2 n QQ 

E E / *ll {H 1 + e-^^W-")] + ln[l + e-^MbnH/*)]] . (27) 



k=0 A=l 
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Integrating by parts with respect to pii in Eq.([Z7|) we find 

£ eff (B, /i, T) = g g ^ dp,, ^-L | _____ + _____ | . 

(28) 

Using l/(l-|-e /3( '' Efe ' A<i 'ii) T ^) — ► #[±/i — -^^(py)], as /3 — > oo, we may in the limit of vanishing 
temperature perform the momentum integral and arrive at Eq . fl24|) . 



with respect to the magnetic field, M = -J^£ c{i , with the result 



The magnetisation of the fermion gas is now easily found by performing the derivative 

_a_ 

dB' 



r_ 

M(B,fi)= _ M n (B,fi) , (29) 

n=0 

where the contribution from the n-th Landau level is 

M n (B, l i) = b n ^ 1 Lyftf - m* - 2eBn - (m 2 + AeBn) In + ^Z"? ~ 5^ ) } . 

(30) 

We notice that lim (At 2_ m 2_ 2e , B?w0 + ) £ n ( J B, \i) = 0, and Yim. tyP 2_ m 2_ 2eBn ^+ ) M n (B = 0, so 
that the Lagrangian density as well as the magnetisation are continuous. Figure 3 shows 
the total magnetisation, M tot (B, fi) = -^C eS (B, ji) ||, as a function of the magnetic field 
for fixed chemical potential (however, the vacuum contribution -Jj^£ cS (B) is small in this 
range of parameters). In agreement with Refs. H, |6], we see that for low temperatures 
(here T = 0), the relativistic fermion gas exhibits the de Haas-van Alphen effect. 
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Figure 1: The fermion density as a function of the chemical potential for fixed magnetic held, 
eB/m 2 = 5. 
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Figure 2: The fermion density as a function of the magnetic held for hxed chemical potential, 
[i = 5. 




Figure 3: The total magnetisation as a function of the magnetic held for hxed chemical 
potential, fj, = 5. 



